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Abstract 



Archimedean copulas are popular in the world o f multivariate modelli ng as a 
result of their breadth, tractability, and flexibility. McNeil Sz Neslehoval (2009) 



showed that the class of Archimedean copulas coincides with the class of multi- 
variate £i-norm symmetric distributions. Building upon their results, we intro- 
duce a class of multivariate Markov processes that we call Archimedean survival 
processes' (ASPs). An ASP is defined over a finite time interval, is equivalent in 
law to a multivariate gamma process, and its terminal value has an Archimedean 
survival copula. There exists a bijection from the class of ASPs to the class 
of Archimedean copulas. We provide various characterisations of ASPs, and a 
generalisation. 

Keywords: Archimedean copula, gamma process, gamma bridge, multivariate 
Liouville distribution 



1 Introduction 

The use of copulas has become commonplace for dependence modelling in finance, 
insurance, and risk management (see, for example, Cherubini et al. jsl]. Frees & Valdez 
j^, and McNeil et al. [l^). The Archimedean copulas — a subclass of copulas — have 
received particular attention in the literature for both their tractability and practical 
convenience. An ra-dimensional Archimedean copula C : [0, 1]" — )■ [0, 1] can be written 
as 

du) = h{h-\u,) + ■ ■ ■ + h-\un)), (1) 

where h is the generator function of C. We introduce a family of multivariate stochastic 
processes that we call Archimedean survival processes (ASPs). ASPs are constructed 
in such a way that they are naturally linked to Archimedean copulas. An ASP is 
defined over a finite time horizon, and its terminal value has an £i-norm symmetric 
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distribution. This implies that the terminal value of an ASP has an Archimedean 
survival copula. Indeed, there is a bijection from the class of Archimedean copulas to 
the class of ASPs. 

A random vector X has a multivariate Liouville distribution if 

^'=R^^, (2) 

where Ris a non-negative random variable, and G is a vector of n independent gamma 
random variables with identical scale parameters (see, for example. Fang et al. [5]). In 
the special case where G is a vector of identical exponential random variables, X has 



an £i-norm symmetric distribution. McNeil & Neslehova [10[ give a account of how 
Archimedean copulas coincide with survival copulas of £i-norm symmetric distributions 
which have no point-mass at the origin. This particular relationship relies on the 
characterization of n-monotone functions through an integral transform of Williamson 



171 ]. Then in iMcNeil Sz Neslehoval generalise Archimedean copulas to so-called 
Liouville copulas, which are defined as the survival copulas of multivariate Liouville 
distributions. 

Norberg [l^ suggested using a randomly-scaled gamma bridge (also called a Dirich- 
let process) for modelling the cumulative payments made on insurance claims (see also 
Brody et al jl]). Such a process {^tT}o<t<T can be constructed as 

^tT = RltT, (3) 

where i? is a positive random variable and {'ytr} is an independent gamma bridge 
satisfying 7or = and 7Tr = 1, for some T G (0, oo). This is an increasing process 
and so lends itself to the modelling of cumulative gains or losses, where the random 
variable R represents the total, final gain. We can int erpret -R as a signal and the 



gamma bridge {jtr} as independent multiplicative noise. iBrody et al\ show that {^tr} 
is a Markov process, and that 

where u is the law of R, and m > is a parameter. The process {^jt} can be considered 
to be a gamma process conditioned to have the marginal law u at time T, and so belongs 
to the class of Levy random bridges (see Hoyle et al. jsj). As such, we call a process 
that can be decomposed as in ([3]) a 'gamma random bridge' (GRB). 

Archimedean survival processes are an ra-dimensional extension of gamma ran- 
dom bridges. Let the process {{^t^\^t^\ ■ ■ ■ ,C,t'^^)~^}o<t<T be an ASP. Then each one- 
dimensional marginal process {^f''} is a GRB. Thus we can write 

#=X.7!?, (5) 

for some gamma bridge {7{'r} and some independent Xi > 0. The Xj's are iden- 
tical but in general not independent, and the {7^7-^ }'s are identical but in general 
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not independent. We shall construct each {^f"*} by splitting a 'master' GRB into 
n non-overlapping subprocesses. This method of splitting a Levy random bridge into 
subprocesses (which are themselves Levy random bridges) was used by Hoyle et al. 
[3] to develop a bivariate insurance reserving model based on random bridges of the 
stable-1/2 subordinator. A remarkable feature of the proposed construction is that the 
terminal vector (i^r^ ■ ■ ■ •> ^t^Y has an £i-norm symmetric distribution, and hence 
an Archimedean survival copula. In particular, we shall show that 



P 

where 



> F{iP) >u„..., > = I V F-\u,) ) , (6) 




(u) = P 



>u , fori = l,2,...,n. (7) 



Here F{x) is the marginal survival function of the ^^^'s, and F~^{u) is its (generalised) 
inverse. The right-hand side of ([6]) is an Archimedean copula with the generator func- 
tion F{x). 

We shall also construct Liouville processes by splitting a GRB into n pieces. By 
allowing more flexibility in the splitting mechanism and by employing some determin- 
istic time changes, a broader range of behaviour can be achieved by Liouville processes 
than ASPs. For example, the one-dimensional marginal processes of Liouville process 
are in general not identical. 

A direct application of ASPs and Liouville processes is to the modelling of multivari- 
ate cumulative gain (or loss) processes. Consider, for example, an insurance company 
that underwrites several lines of motor business (such as personal motor, fleet motor 
or private-hire vehicles) for a given accident year. A substantial payment made on one 
line of business is unlikely to coincide with a substantial payment made on another line 
of business (e.g. a large payment is unlikely to be made on a personal motor claim at 
the same time as a large payment is made on a fleet motor claim). However, the total 
sums of claims arising from the lines of business will depend on certain common factors 
such as prolonged periods of adverse weather or the quality of the underwriting process 
at the company. Such common factors will produce dependence across the lines. An 
ASP or a Liouville process might be a suitable model for the cumulative paid-claims 
processes of the lines of motor business. The one-dimensional marginal processes of a 
Liouville process are increasing and do not exhibit simultaneous large jumps, but they 
can display strong correlation. 

ASPs can be used to interpolate the dependence structure when using Archimedean 
copulas in discrete-time models. Consider a risk model where the marginal distributions 
of the returns on n assets are fitted for the future dates ti < ■ ■ ■ < t„ < T < oo. An 
Archimedean copula C is used to model the dependence of the returns to time T. At 
this stage we have a model for the joint distribution of returns to time T, but we have 
only the one-dimensional marginal distributions at the intertemporal times ti, . . . ,tn- 
The problem then is to choose copulas to complete the joint distributions of the returns 
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to the times ti, . . . in a way that is consistent with the time-T joint distribution. 
For each time ti, this can achieved by using the the time-tj survival copula implied by 
the Archimedean survival process with survival copula C at terminal time T. 

This paper is organized as follows: In Section 2, we review multivariate £i-norm 
symmetric distributions, multivariate Liouville distributions, Archimedean copulas and 
gamma random bridges. In Section 3, we define ASPs and provide various characteri- 
sations of their law. We also detail how to construct a multivariate process such that 
each one- dimensional marginal is uniformly distribution. In Section 4, we generalise 
ASPs to Liouville processes. 



2 Preliminaries 

This work draws together ideas from mathematical statistics and the theory of stochas- 
tic processes. This extended preliminary section gives relevant background results from 
both of these subjects. 

We fix a probability space {Q, P, J-") and assume that all processes and filtrations un- 
der consideration are cadlag. We let denote the generalised inverse of a monotonic 
function /, i.e. 

, , I infix : f(x) > y}, f increasing, , ^ 

r\y) = { . ,] \)[- , ' . ^' (8) 

I mtjx j{x) <y\, J decreasmg. 
We denote the li norm of a vector x G M" by ||x||, i.e. 

n 

M=Y.\x,\. (9) 



2.1 Multivariate distributions 

In this subsection we present some definitions and results from the theory of multivari- 
ate distributions. We refer the reader to the thorough exposition by Fang et al. ^ for 
further details. 

2.1.1 Multivariate £i-norm symmetric distributions 

The multivariate £i-norm symmetric distributions form a family of distributions that 
are closely related to Archimedean copulas. The definition of the n-dimensional ^l- 
norm symmetric distribution is in terms of a random variable uniformly distributed on 
the simplex 

5 = {uG [0,1]" : ||u|| =1}. (10) 
Such a random variable U has the stochastic representation 
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where E is a vector of n independent, identically-distributed, exponential random vari- 
ables. Note that this representation holds for any value of the rate parameter A > 
of the exponential random variables, and that the random variable ||E|| has a gamma 
distribution with shape parameter n, and scale parameter A~^. Each marginal variable 
Ui has a beta distribution with parameters a = 1 and /3 = n — 1; thus the survival 
function of Ui is 

F[Ui>u] = {l-ur-\ (12) 

for < M < 1. 

Definition 2.1. A random variable X taking values in M" has a multivariate £i-norm 
symmetric distribution if 

X '= KU, (13) 

where R is a non-negative random variable, and \J is a random vector uniformly dis- 
tributed on the simplex S. We say that the law of R is the generating law of the 
distribution. 

Remark 2.2. The construction of multivariate ii-norm symmetric random variables 
is similar to the construction of elliptical random variables. To be precise, in f lT3|) if 
U was uniformly distributed on the unit sphere in M", then X would have an elliptical 
distribution. 

Note that if R admits a density, then X satisfying ( 1T3|) admits a density, and this 
density is simplectly contoured. This is analogous to the elliptical contours of elliptical 
distributions. 

If X is a multivariate £i-norm symmetric random variable with generating law v, 
then the survival function of each one- dimensional marginal of X is 

F{x) = r[Xi > x] 

'J 

(1 -x/r)"-V(dr), (14) 

for X > 0. The survival function F determines the law u. Indeed, using the results of 
Williamson 17 , McNeil & Neslehova 10 showed that 



^(-IjWri^ (-lj--x--max[U,i^o- ^'{x)\ 
^^([U,xjj-1 2^ ^, , (15) 

where F^^^ is the fcth derivative of F, and 

Fo X = <^ -, . 16 
^ ^ [ 1 -F(0), X = 0. ^ ' 

The following theorem provides the multivariate version of f|T4l) : a proof can be found 
in Fang et al. [l, Theorem 5.4]. 
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Theorem 2.3. //X has a multivariate ii-norm symmetric distribution with generating 
law V, then the joint survival function of~K. is 



/oo 
(1- ||x||/rr-V(dr) 
..xll 



Ffllxl 
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/or X e Wl. 



2.1.2 Multivariate Liouville distributions 

The muhivariate Liouville distribution is an extension of the multivariate ^i-norm sym- 
metric distribution. Before defining the multivariate Liouville distribution, it is conve- 
nient to first define the Dirichlet distribution. The n-dimensional Dirichlet distribution 
is a distribution on the simplex S defined in ffTOjl . 



Definition 2.4. Let G be vector of independent random variables such that Gi is a 
gamma random variable with shape parameter ai > and scale parameter unity. Then 
the random vector 

has a Dirichlet distribution with parameter vector o: = (ai, . . . , On)"''- 

Remark 2.5. The scaling property of the gamma distribution implies that kG, k > 0, 
is a vector of gamma random variables each with scale parameter k. Since (|T7I) holds, 
if we replace G with nG, we could have used an arbitrary positive scale parameter in 



Definition \2.4 



In two dimensions, a Dirichlet random variable can be written as [B, 1 — B)~^ , 
where S is a beta random variable. If all the elements of the parameter vector a 
are identical, then D is said to have a symmetric Dirichlet distribution. Notice that 
if ctj = 1 for i = 1,2, ... ,n, then D is uniformly distributed on the simplex S. The 
density of {Di, D2, . . . , Dn-i)~^ is 



X 1-^ 



for X G [0, 1]" \ ||x|| < 1, where Xn = 1 — Yl^=i T[z] is the gamma function, 

defined as usual for x > by 



00 

r[a;] = / M^-^e-"dM. (19) 
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The first- and second-order moments of the Dirichlet distribution are given by 

rv ■ 

E[A] = (20) 
ll^ll 

Varlfl.l = ,°''!"^:?-' . (21) 

||Q;||^(||a|| + 1) 

°°-l^'-^''-- ||c|P(IM + i) - 

The Dirichlet distribution is an extension of a random variable uniformly distributed 
on a simplex. The multivariate Liouville distribution is a similar extension of the 
multivariate £i-norm symmetric distribution. 

Definition 2.6. A random variable X has a multivariate Liouville distribution if 

X ^= i?D, (23) 

for R > a random variable, and D a Dirichlet random variable with parameter 
vector a. We call the law of R the generating law and a the parameter vector of the 
distribution. 



^11 iiJni„iiyi"ll-i ii r[ai] ^ ^ 



In the case where R has a density p, the density of X exists and can be written as 

,iNllryll-1 I 

for X G M". Writing /xi = K[R] and /i2 = IE[i?^] (when these moments exist), the first- 
and second-order moments of X are given by 

ry ■ 

E[X,]=/ii^, (25) 
||q:|| 

\\a\\ \ CK + i CK / 



Cov[X., X,] = ^ ( - y for . ^ (27) 

||q;|| \ ||q;|| + i ||q;|| / 

2.2 Archimedean copulas 

A copula is a distribution function on the unit hypercube with the added property that 
each one-dimensional marginal distribution is uniform. For further details, we refer to 



Nelsen [ij]. We define a copula as follows: 



Definition 2.7. An n-copula defined on the n-dimensional unit hypercube [0, 1]" is a 
function C : [0, 1]" — )■ [0, 1], which satisfies the following: 

1. C{u) = whenever uj = for at least one j = 1, 2, .., n. 
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2. C{u) = Uj if Ui = 1 for all i ^ j. 

3. C is n-increasing on [0, 1]", that is 
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5^ ■ ■ ■ . . . , «„,J > 0, (28) 

for all U2,i, . . . , Un,iy and (■Ui,2, ^2,2, • • • , Un,2V in [0, 1]" with uj^i < Uj^2- 

In the definition above, condition [3] is necessary to ensure that the function C is a 
well-defined distribution function. The theory of copulas is founded upon a theorem 
of Sklar. This theorem was reformulated in terms of survival functions by McNeil & 
Neslehova 10| as follows: 

Theorem 2.8. Let H he an n- dimensional survival function with margins Fi, i = 
1,2, ... ,n. Then there exists a copula C , referred to as the survival copula of H , such 
that, for any x G M", 

^(x)=C(Fi(xi),...,F„(x„)). (29) 
Furthermore, C is uniquely determined on 

D = [0, 1]" : M G ranFi x ■ ■ ■ x ranF„} , 

where ran/ denotes the range of f . In addition, for any u E D, 

C{u) = H{F{\u,),...,F~\un)). 

Conversely, given a copula C and univariate survival functions Fi, i = 1, . . . ,n, H 
defined by fl2^ is an n-dimensional survival function with marginals Fi, . . . ,Fn and 
survival copula C . 

From a modelling perspective, one of the attractive features of copulas is that they 
allow the fitting of one- dimensional marginal distributions to be performed separately 
from the fitting of cross-sectional dependence. However, this two-step approach of 
modelling multivariate phenomena by first specifying marginals and then choosing a 
copula is not suited to all situations (for criticism see, for example, Mikosch [lit). 

Archimedean copulas are copulas that take a particular functional form. The fol- 
lowing definition given in 10|] is convenient for the present work: 

Definition 2.9. A decreasing and continuous function h : [0, oo) — )■ [0, 1] which sat- 
isfies the conditions h{0) = 1 and lima;^oo h{x) = 0, and is strictly decreasing on 
[0,inf{x : h{x) = 0}] is called an Archimedean generator. An n-dimensional copula C 
is called an Archimedean copula if it permits the representation 

Civi) = h{h-\u,) + ■■■ + h-\Un)), U G [0, 1]", 

for some Archimedean generator h with inverse h"^ : [0,1] — i- [0,oo), where we set 
h{oo) = and h-\0) = ml{u : h{u) = 0}. 
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If X is a random vector with a multivariate ^i-norm symmetric distribution such 
that P[X = 0] = 0, then its marginal survival function F given in f[T^ is continuous. 
Hence it follows from Theorem 12.31 that 



P[F(Xi) > m, F{X2) >U2,..., F{Xn) > n„] = F V F-\u,) . (30) 




In other words, X has an Archimedean survival copula with generating function h{x) = 
F{x). McNeil & Neslehova 10| showed that the converse is also true: 

Theorem 2.10. Let U foe a random vector whose distribution function is an n-dimensional 
Archimedean copula C with generator h. Then {h~^ (Ui) , h^^ {U2) , ■ ■ ■ , h~^ {Un))~^ has 
a multivariate ti-norm distribution with survival copula C and generating law v . Fur- 
thermore, V is uniquely determined by 

^ (-l)"-ix"-^max[0,/i("-^)(x)] 

fc=0 ^ ^ 

Remark 2.11. There is one-to-one mapping from distribution functions on the posi- 
tive half-line to the class of n-dimensional Archimedean copulas through the invertible 
transformation u ^ h. 



2.3 Gamma random bridges 

A gamma random bridge is an increasing stochastic process, and both the gamma 
process and gamma bridge are special cases. 



2.3.1 Gamma process 

A gamma process is a subordinator (an increasing Levy process) with gamma dis- 
tributed increments (see, for example, Sato [l6[). The law of a gamma process is 
uniquely determined by its mean and variance at time 1, which are both positive. Let 
{'jt} be a gamma process with mean and variance m > at time 1; then 

E[7t] = mt, and Var[7j] = mt. (31) 

The density of ■jt is 

m = l{-o}^e-^ (32) 

Due to the scaling property of the gamma distribution, if k > then the process {K'jt} 
is a gamma process with mean mn, and variance m/t^ at t = 1. The characteristic 
function of 74 is 

E[e'^^*] = (1 - iA)^™*, for A G C. (33) 
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As noted in Brody et al. [2.], the parameter m has units of inverse time, and so {7f} 
is dimensionless. Taking k = 1/m, the scaled process {^7*} has units of time, making 
this ahernative parameterisation suitable as a basis for a stochastic time change (see, 
for example, Madan & Seneta joj). The characteristic function of K'yt is then 

Efe^^'^^'] = (1 -iA/m)-"^*. (34) 

In the limit m — )■ oo this characteristic function is e''^*, which is the characteristic 
function of the Dirac measure centred at t. It follows that {^7^} {t} as m — )■ 00. 



2.3.2 Gamma bridge 

A gamma bridge is a gamma process conditioned to have a fixed value at a fixed future 
time. A gamma bridge is a Levy bridge, and hence a Markov process. Gamma bridges 
exhibit a number of remarkable similarities to Brownian bridges, some of which have 
been presented by Emery & Yor [3]. Let {'ytT}o<t<T be a gamma bridge identical in law 
to the gamma process {jt} pinned to the value 1 at time T. Using the Bayes theorem, 
the transition law of {'jtr} is given by 

P [jtT e dyljsT = x] =F [7i G dy I 7^ = x, 7t = 1] 

ft-s{y - x)fT-t{l -y) 



/t-.(I-x) 

{l-x)B[m{t-s),m{T - t)] 



for < s < t < T and a; > 0. Here B[a, /3] is the beta function, defined for a, /3 > by 
B[a,P]= [\-\l-xr'dx = ^jf^^. (36) 

Jo r[« + p\ 

We say that m is the activity parameter of {'jtr}- If the gamma bridge {jtr} has reached 
the value x at time s, then it must yet travel a distance 1 — x over the time period 
(s, T]. Equation (135|) shows that the proportion of this distance that the gamma bridge 
will cover over {s,t] is a random variable with a beta distribution (with parameters 
a = m{t — s) and f3 = m{T — t)) . The conditional characteristic function of 'jtr is 

E [e^^^*^ I 7,T = = M[m{t - s),m{T - s),i{l - x)X\, (37) 

where M[a, (3, z] is Kummer's confluent hypergeometric function of the first kind, which 
can be expanded as the power series [H, 13.1.2] 

, ,r ^ 1 a a(a + 1) a(a + l)(a + 2) , . 

M[a./;..] = l + -. + A__|_+ ^|^^^n^^^> _ + .... (38) 
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Taking the hmit as m — )■ oo in fl37|) . we have 



k=0 



t-sV (i(l-x)A)' 



T-s k\ 



= exp(^i^-^(l-x)A^ , (39) 

which is the characteristic function of the Dirac measure centered at {1 — x){t — s) / [T — 

s). It then follows from the Markov property of gamma bridges that {^tr} ^ {t/T} 
as m — )■ oo. 

It is a property of gamma processes that the renormalised process {'jt/lT}o<t<T 
is independent of jt (indeed, this independence property characterises the gamma 
process among Levy processes). This leads to the remarkable identity 

^] = Ht}. (40) 

The identity f HOj) can be proved by showing that the process on the left-hand side is 
Markov, and then verifying that its transition law is the same as f l35|) . This can be done 
using the results in Brody et al. [2]. We note three properties of gamma bridges that 
follow from fHU]) . The first is that the bridge of the scaled gamma process {^7^} is, for 
any k > 0, identical in law to the bridge of the unsealed process. The second property 
is that a {7t}-bridge to the value 2; > at time T is identical in law to the process 
{z'^it}- The third is that the joint distribution of increments of a gamma bridge is 
Dirichlet. To see this last fact, fix times = to < < ■ ■ ■ < = ^ and define 

^^=lU-lU_,, (41) 

Ai = -iu,T - lu^uT- (42) 

Then Aj has a gamma distribution with shape parameter = m{ti — ti_i) and scale 
parameter unity. Hence 

A A ^ (Ai,A2,...,A^) 

The fact then follows from Definition 12.41 



2.3.3 Gamma random bridge 

" ' as 



Under a different name, the gamma random bridge was introduced by Norberg [15 
a model for a cumulative payment process in an insurance model. It was also studied 
in detail by Brody et al. [2!] . We define a gamma random bridge as follows: 

Definition 2.12. The process {Tt}o<t<T is a gamma random bridge if 

{T,} {R^,t}, (44) 
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for R > a random variable, and {'ytr} o. gamma bridge. We say that {Tf} has 
generating law u and activity parameter m, where of v is the law of R and m is the 
activity parameter of {'ytr}- 

Remark 2.13. Suppose that {Tt} is a GRB satisfying (HH). // P[i? = z] = I for 
some z > {], then {Tt} is a gamma bridge. If R is gamma random variable with 
shape parameter niT and scale parameter n, then {Tt} is a gamma process such that 
E[Tt] = ruKt and Yax\Vt] = tukH, for t G [0,T). 

Gamma random bridges (GRBs) fall within the class of Levy random bridges de- 
scribed by Hoyle et al. [Sj. The process {r^} is identical in law to a gamma process 
defined over [0, T] conditioned to have the law of R at time T. The bridges of a GRB 
are gamma bridges. GRBs are Markov processes, and the transition law of {Ff} is 
given by 



F[Tt edy\Ts = x 

^ hi 

B[m(T - t),m{t - s)] J^{z - x)'^'(T-s)-i^i-mT j^^^z 



and 

p|r. e d, I r. ^ .] = '£> , (46) 

where B[a, /3] is the beta function. Since increments of a gamma bridge have a Dirichlet 
distribution, it follows from Definition 12.61 that the increments of a gamma random 
bridge have a multivariate Liouville distribution. The following proposition, stated as 
a corollary in jsj for a general Levy random bridge, is a key result for the construction 
of Archimedean survival processes: 

Proposition 2.14. Let {F^} be a GRB with terminal law v and activity parameter m. 

(A) Fix times si,Ti satisfying < Ti < T — si. The time-shifted, space- shifted partial 
process 

d'^ = F,,+i-F,„ (0<t<Ti), 
is a gamma random bridge with activity parameter m, and with generating law 

mT-i-l roo 

i/«(da;) = / z'^^-\z - xT^''-^'^-^p{dz)dx. 



B[mTi,m(T-Ti)] 

(B) Construct partial processes {^t^}Q<t<T^, i = 1, . . . ,n, from non- overlapping por- 
tions of {Tt}, in a similar way to that above. The intervals [si,Si + Ti], i = 
l,...,n, are non- overlapping except possibly at the endpoints. Set = ^^.^ 
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when t > Ti. If u> t, then 



P 



(n) _ An) 



P 



ei')-d'^<xi,...,er-er<x. 



in) _ f (n) 



(0 



?■ = 1 2 



where the filtration {Tt\ is is given by 
Remark 2.15. Define the process {Rt} by 

n 



(47) 



i=l 



for t G [0,maXjTj]. Then {Rt} is a GRB with terminal law v, and time- dependent 
activity parameter 



M{t)=mY,Mt<T.}- 



(4J 



j=i 



The proof of this result is similar to the proof of the special case that appears later in 
Proposition 13.61 

We can construct an ra-dimensional Markov process {$,^} from the partial processes 
of Lemma B. 141 part (B), by setting 



~ ) • • • ) St ) 



(")^T 



(49) 



The Markov property means that, for any fixed time s > 0, the J-'^-conditional law 
of {^t}s<t is identical to the ^^-conditional law of {^t}s<t. The remarkable feature of 
Lemma 12.141 part (B), is that the J-'^-conditional law of {^^ — ^s}s<t is identical to 
the i?<j-conditional law of {^^ — ^s}s<t. Hence the increment probabilities of the n- 
dimensional process {^^} can be described by the one-dimensional state process {Rt}- 



3 Archimedean survival process 

We construct an Archimedean survival process (ASP) by splitting a gamma random 
bridge into n non-overlapping subprocesses. We start with a 'master' GRB {Tt}Q<t<n 
with activity parameter m = 1 and generating law z/, where n G N+, n > 2. In 
this section, we write ft for the gamma density with shape parameter unity and scale 
parameter unity (in fl32|) we set m = 1). That is 

m = -Ywr- (50) 
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Definition 3.1. The process {^t}o<t<i "^s an n-dimensional Archimedean survival pro- 
cess if 



r ^ -1 law 
l€tJO<t<l = < 



\_^(n-l)+t 



n-1 



0<t<l 



for {rt}o<t<n is a gamma random bridge with activity parameter m = 1. We say that 
the generating law of {Tt} is the generating law o/j^^}. 



Note that from Definition 12121 P[r„ = 0] = 0, and so P[|i = 0] = 0. Each one- 
dimensional marginal process of an ASP is a subprocess of a GRB, and hence a GRB. 
Thus ASPs are a multivariate generalisation of GRBs. We defined ASPs over the time 
interval [0, 1]; it is straightforward to restate the definition to cover an arbitrary closed 
interval. 

Proposition 3.2. The terminal value of an Archimedean survival process has an 
Archimedean survival copula. 

Proof. Let {$,^} be an n-dimensional ASP with generating law u. Then we have 
P[^i G dx] = P [Fi G dxi, - Fi G dxs, . . . , F„ - F„„i G dx„] 
= P 



K — G dxi, K G ax2, 



T-)1n 1n~l ^ 1 

H G axr 



In 



In 



In 



(51) 



for X G M", R a random variable with law z/, and {7f} a gamma process such that 7f 
has the density (l50l) . Each increment 7j — 7j_i has an exponential distribution (with 
unit rate). Thus 



P[^i G dx] = P 



E 

G dx 

E 



(52) 



for E an n- vector of independent, identically-distributed, exponential random variables. 
From Definition 12.11 has a multivariate £i-norm symmetric distribution. Therefore, 
it has an Archimedean survival copula. 



Remark 3.3. Let gi : ]R_|_ — )■ M 6e strictly decreasing for i 
an ASP. Then the vector-valued process 



□ 



1, . . . ,n, and let {^^j be 



^5i(d'^),...,^.(rf^),...,^n(d"^; 

has an Archimedean copula at time t = 1. 



'} 



0<t<l 
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3.1 Characterisations 

In this subsection we shall characterize ASPs first through their finite-dimensional 
distributions, and then through their transition probabilities. 

3.1.1 Finite-dimensional distributions 

The finite-dimensional distributions of the master process {r^} are given by 

F[Tt^ G dxi, . . . , G dxk, r„ G d^] = F[Tt^ G dxi, . . . , T^^ G dxk | r„ = u{dz), 

(53) 

where Xq = 0, for all k G N+, all partitions = to < < ■ ■ • < ^fc < all 2; G M+, 
and all (xi, . . . ,Xk)~^ = x G M^. It was mentioned earlier that the bridges of a GRB 
are gamma bridges. (In fact, this is the basis of the definition of Levy random bridges 
given in Hoyle et al. [8|.) Hence, for {7^} a gamma process such that ]E[7i] = 1 and 
Var[7i] = 1, we have 

W{Tt, G dxi,...,Vt, G dxfc,r„ G dz] 

= P[7t, G dxi, . . . , 7f^ G dxk \jn = z] i^{dz). (54) 

From (HOD and dSj), we have 

1 R 

{Tt, - Tto,...,Tt, - Tt,_,,Tn-TtJ = — (7ii -7io.---.7ife - 7ife-i, 7n - 7iJ- (55) 

In 

Hence, from Definition 12. 6^ (F^^ — F^q, . . . , F^^. — Tt^_^, F„ — F^^)^ has a multivariate 
Liouville distribution with generating law u and parameter vector (ti — to, . . . ,tk — 

tk-i, n — tkY . 

We can use these results to characterise the law of the ASP {^^} through the joint 
distribution of its increments. Fix ki > 1 and the partitions 

= 4<'^1<---<4, = 1, (56) 
for i = 1, . . . ,n. Then define the non-overlapping increments {Ajj} by 

\j = - , (57) 
for j = 1, . . . , ki and i = 1, . . . ,n. The distribution of the vector 

A = (An, A12,..., Aifc,, 

A2I, A22, • • • , A2A;2, 

A„i, A„2, • • • , A^fc^)"'' (58) 

characterises the finite-dimensional distributions of the ASP {$t}- Thus it follows from 
the Kolmogorov extension theorem that the distribution of A characterises the law of 
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{^t}. Note that A contains non-overlapping increments of the master GRB {F^} such 
that ||A|| = r^. Hence A has a multivariate Liouville distribution with parameter 
vector 



J.2 J.2 J.2 J.2 J.2 J.2 

''1 ~ ''0) '-2 ~ ''15 • • • 5 ''A;2 ~ ''A;2-l' 



^1 '^0?^2 ^l5---)^fc„ ^L-l)"'"' (^^) 



and the generating law i/. 



3.1.2 Transition law 

We denote the filtration generated by {^t}o<t<i by {J^t}- From Lemma l!^. 141 {^j} is a 
Markov process with respect to {J-'t}- We shall calculate the transition probabilities of 
{^f-} after introducing some further notation. 

For a set 5 C M and a constant x G M, we write B + x for the shifted set 

B + x = {y eR:y -X e B}. (60) 

In what follows, we assume that {^^j is an n-dimensional ASP with generating law u, 
and that {F^} is a master process of {$,t}- We define the process {Rt}o<t<i by setting 

n 

^* = E^*'^ = 11^*11- (61) 

i=l 

Note that the terminal value of {Rt} is the terminal value of the master process {F^}, 
i.e. Ri = F„. We define a family of unnormalised measures, indexed by t G [0, 1) and 
X G M+, as follows: 

MB;x) = uiB), (62) 

Jb fn{z) 



F[n]e^ 



l{.>.}^'-"(^-x)"(i-*)-iz/(dz) (63) 



r[^(i-t)] Jb 

for 5 G B{M.). We also write 

^i(a;) = V't([0,oo);x). (64) 
It follows from f l5^ and the independent increments of gamma processes that 

F[rt, G dxi, . . . , Ft, G dxk, Tn G dz] = l[[f u-u_,{x^ - x,_i) dx,] ^""^^.^^ ~ ^"^ //(dz) 

fc 

= n['^*'-*>-i(^* -a;i_i)dxi]^/'t,/„(d2;;a;„). (65) 



1=1 
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Proposition 3.4. The ASP {^,^} is a Markov process with the transition law given by 



P 



.(n-l) 



(n) 



^.fllxl 



r[i-s] 



and 



i=l 



(67) 



where T{t) = I - {I -t)/n, < s <t < I, and B e i3(M). 

Proof. We begin by verifying flBBl) . From the Bayes theorem we have 



P 



p 



G d^i, . . . , er G d^„_i, ll^ill G 5 + L e d 



.(n-1) 



n-1 



P G dx] 

The law of Ri = is u; hence using (165|) the numerator of (l68l) is 



(6^ 



P 



ef^Gdzi,...,er'^Gd^„_i,CG dx 



Ri = u 



n-l 



n[/s(x.)dx.]n[/i-.( 

1=1 



and the denominator is 

n 

P e d^\Ri=u] U{du) = lllfsiXi) dXi] 



u=0 



1=1 



fi-sju - Er=i 

/n(l-.)(^^ - ||x||) 



u=0 



fniu) 



u{du), (69) 



z/(du). (70) 



In equations fl69|) and flTOj) we have used the fact that, given = Ri, {$,t} is a vector 
of subprocesses of a gamma bridge. Dividing f l69|) by flTOj) yields 



r^D,v-n-i -T-TT fi-sf^^ — y^r 1 -Zj) t^(d?i) ^_^} 



L=o - II^ID'^W 



i=l 



i^ris){B + 2;^; x„ + j-j (^»-x^) "e ^^^^ 



^s([0,oo); ||x| 



i=l 



T[l-s] 



as required. 

We shall now verify (E7|) following similar steps. From the Bayes theorem we have 



P[^, G dy I = x] 



P[^, Gdy.^Gdx] 
ms e dx] 



(72) 
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The numerator of (1721) is 



/•oo 

/ P[^tGdy,CGdx|i?i = z] v{dz) = 

Jz=0 

flifsix,) dx,] flUUV: - X,) dy,] r /-a-^)(^ -||yll) ^(dz), (73) 
i=l i=l -^^=0 •^'^^^^ 

and the denominator is given in (!70l) . Dividing ( !73|) by (ITOll yields 

TOO '^"i ^ . 7 []— [7T — rrr TT[/t-s(l/i ~ ^j) dy?,] = 

V^i([0, oo); lly|l) A (y. - a:,)ft--)-ie-(^'-"') 
^.([O,oo);||x||)ll r[t-s] ^'^^ 



which completes the proof. □ 
Remark 3.5. When the generating law v admits a density p, (168|) is equivalent to 



3.1.3 Increments of ASPs 

We shall show that increments of an ASP have n-dimensional Liouville distributions. 
Indeed, at time s G [0, 1), the increment — t G (s, 1], has a multivariate Liouville 
distribution with a generating law that can be expressed in terms of the ^^-conditional 
law of the norm variable Rt = ||CJ- Before we show this, we first examine the law of 
the process {Rt}- 

Proposition 3.6. The process \Rt}o<t<T is a GRB with generating law v and activity 
parameter n. That is, 

F[Rt G dr I 

and 

F[Ri edr\$^ 
forO<s<t<l. 

Before proceeding the proof, note that, after simplification, (176|) and (177|) are con- 
sistent with (gSD and (06]). 



^ vl>,(r) (r- ||x||)"(*--)-^exp(-(r- ||x||)) ^^ ^^^^ 



X 



T[n{t - s)] 



^s(dr; ||x| 
*s(||x||) 



(77) 



18 



Proof. Since {^j} is a Markov process with respect to {J^t}, {Rt} is a Markov process 
with respect to {J^t}- Thus to prove the proposition we need only verify that the 
transition probabihties of {Rt} match those given in fl7^ and fITTI) . 

We first verify the ^^-conditional law of Ri. We can calculate this using the Bayes 
theorem; 



_ ;^/n(i-s)(r- ||x||)^(dr) 
V^s(dr; ||x||) 

Similarly, the ^^-conditional law of Rt for t G (s, 1) can be derived by the use of 
the Bayes theorem; 

JZo e dx, Rtedr\R, = z]F[R, e d^] 



¥[Rt e dr I = x] 



JZo IZo ^i^s e dx, Rtedr\Ri = z]drF[R^ G d^] 
IZo ■J^)fn{t-s){r - ||x||)/„(i_t)(2; - r)dr p{dz) 
IZo A^fe) Xl||x|| fn{t-s){r - M)fn(i-t){z - r)dr iy{dz) 



(79) 



-/„(t_,)(r-||x||)dr. (80) 



The denominator of (1791) was simplified using the fact that gamma densities with 
common scale parameter are closed under convolution. □ 

We define the measure i^st, < s < t < 1, by 

iyst{B) = nRteB\Q. (81) 

Thus we have 

^ ips{dr]Rs) 

and 

^ ~ ^JrZ) nn{t-s)] < ^^^> 

When Ust admits a density, we denote it by Pst{r) = i'st{dr)/dr. We see from fl83l) that 
Pst exists for t < 1. When t = 1, it follows from the definition of tpt that pst only exists 
if admits a density. 

Note that f[Rt G dr | = F[Rt G dr | R^] for t G (s, 1]. This is not surprising since 
{Rs} is a GRB, and hence a Markov process with respect to its natural filtration. 

We now show that the increments of an ASP are of multivariate Liouville-type. 
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Proposition 3.7. Fix s G [0, 1). Given the increment — t G (s, 1], has an 

n-variate Liouville distribution with generating law 



iy*{B) = u,t{B + R,), (84) 
and parameter vector a = [t — s, . . . ,t — s)^ . 



Proof. First we prove the case t < 1. In this case the density Pst exists. From (1671) and 
we have 

-Ledy I a = 11 T[t - s] 

||y||n(t-.)-l 11 Y[t-s] ^ ' 

Comparing fl85|l to f l2^ shows it to be the law of Liouville distribution with generating 
law Pst{.x + Rs)dx and parameter vector (t — s, . . . ,t — sY . Noting that Pst{.x + Rs)dx = 
u*{dx), where u* is given by flMl) . yields the required result. 

We now consider the case t = 1 when u admits a density p. In this case the density 
Psi exists. From ( I75|) and ( 182|) . we have 



. ^ , T[n]e^^p{\\y\\+Rs) 

[€i G dy I - ^^(^^)(||y|| + 11 r[i _ .1 



s 

T[n{l-t)]p,^i\\y\\+R,)J\ y: 



IL g 

rinrr^d,. (86) 



Hence ^^ — has the required density. 

For the final case where t = 1 and v has no density we only outline the proof since 
the details are far from illuminating. Given the law of — is characterised 
by fl66p . We then need to show that this law is equal to the law of XD, where X is 
a random variable with law v* given by fl84p . and D is a Dirichlet random variable, 
independent of X, with parameter vector (1 — s, . . . , 1 — s)^. This is possible by mixing 
the Dirichlet density f|T8l) with the random scale parameter X. 

□ 

3.2 Moments 

In this subsection we fix a time s G [0, 1), and we assume that the first two moments 
of u exist and are finite. 
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Proposition 3.8. The first- and second-order moments of t G (s, 1], are 



1. 
2. 

3. 



E 



(0 



n 



Var 
Gov 



St 1 St 



1 

n 



t- s + l 



n{t - s) + 1 



fi2 /ii 

n 



t - s 



n 



1^2 



n{t — s) + 1 n{t — s) 



where 



/ii 

At2 



^— ^ {E[Ri I R,] - Rs 

(t-g)(l + n(t-g)) 
(l-s)(l + n(l-s)) 



E[(i?i-i?,)2|i?,]. 



Proof. Given the increment ^t^^s has an n-dimensional Liouville distribution with 
generating law 

u*{A) = Ust{A + Rs), (87) 
for t G (s, 1], and witli parameter vector (t — s, . . . ,t — s)~^ . We have 



poo poo 

/ii = / y u*{dy) = / y ustidy) - i?, = E[i?t 1 1.] - Rs 

Jo J Rs 



(88) 



and 



/i2 



r^*(dy) 



(2/-i?s)'^^.t(dy) = E[(i?t-i?.)'|C 



^9) 



Rs 



It then follows from equations (l25|) - (l27|) that 



1. E 



(i) 



(E[i?,iij-i?,) + ei*), 



2. Var 

3. Gov 



Ai) Aj) 

St 5 St 



1 

n 



t-s + 1 



[E[{R,-R,f\^^]) {E[R,\U-R, 



n{t - s) + 1 



n{t — s) 



To compute E[Rt \ ^g] and E[(i?t — i?^)^ | ^g], we use two results results about Levy 
random bridges found in [sj. First, we can write 

E[Rt I Rs] = I Rs] + ]^Rs- 
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The expression for then follows directly. Second, given Rg, the process {Rt—Rs}s<t<i 
is a GRB with terminal law = i'si{B + Rg) and activity parameter n. Hence, 

given Rg, 

{Rt - Rs}s<t<i = {Xjn}s<t<i, (90) 

where X is a random variable with law z/, and {■yti}s<t<i is a gamma bridge with 
activity parameter n, independent of X, satisfying 7^1 = and 711 = 1. Note that 7^1, 
t G (s, 1), is a beta random variable with parameters a = n{t — s) and /3 = n(l — t). 
Thus 



E[(i?i - Rgf I Rg] = £[72 ]E[X^] 

/■oo 

= E[7,^,] / x^P(dx) 

JO 



/•oo 

E[7'i] / (z/-i?.)'i^.i(dx) 

jRs 

^{R^-Rgf\Rg]. (91) 



(t-5)(l + n(t-s))^.,^ ^ ,2 



(l-s)(l + n(l-s)) 

The expression for 112 follows. □ 
3.3 Measure change 

In this section we shall show that the law of an n-dimensional ASP is equivalent to 
an n-dimensional gamma process. To demonstrate this result it is convenient to begin 
by assuming that under some measure Q the process {$,^} is an n-dimensional gamma 
process, and then show that is an ASP under an equivalent measure P. 
Under Q, we assume that {^^} is an n-dimensional gamma process such that 

n i_i 

Q[ltedx] = n|ij^e--'da;,. (92) 

Hence the gamma processes {^f i = 1, 2, . . . , n, are independent and identical in law. 
The process {-Rt}o<t<i, defined as above by Rt — is a one-dimensional gamma 

process and satisfies 

j.nt—1 

Q[i?tedx] = ^e-^dx. (93) 

As before, the filtration {Tt\ is that generated by 

We shall show that the process {^t(-Rt)}o<t<i is a martingale, where 

,1, iu\ fn{i-t){z - Rt) 
^t[Rt)^ TT-\ ^(d^) 

jRt Jn\Z) 
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For times 0<s<t<l, we have 



Eq [%{Rt)\J^s]=^i 

= Er 



oo 



fn(l-t){.Z - Rt) 



°° fn(l-t){.Z — -Rs — {Rt — Rs)) 



Rt 

oo /"OO 



fn{z) 

fn{i-t)iz - Rs-y) 



u{dz) 



y=0 J z=Rs+y 
oo 



fn{z) 



l^{(iz) fn{t-s){y) 



-Rs 



z=Rs fn{z) Jy=o 
°° fn(l-s)iz — R, 



Rs 



fn{z) 



fn{i-t){z - Rs- y)fn{t-s){y) <lyv{<lz) 

^ v{dz) 



■^s{Rs)- 



(95) 



Since \l/o(-Ro) = 1 and ^t{,Rt) > 0, the process {'^t{Rt)}o<t<i is a Radon-Nikodym 
density process. 



Proposition 3.9. Define a measure P by 

dP 



^^tiRt). 



(96) 



Under P, {^(}o<t<i ^-5 o-'n ASP with generating law v. 



Proof. We prove the proposition by verifying that the transition law of {^^j under P 
is that of an ASP. 

P[|, Gdx|J-,]=Ep[l{C, edx}|J-,] 

= ^^EQ[vI/,(i?,)l{^,Gdx}||,] 

^.(i?t) fr (x.-d^V^-)-V(-'-^"\ 

~^s{Rs)\}, T[t-s] 

Comparing equations ( 1971) and flB7|) completes the proof. 

We can restate the results of this subsection as the following: 



(97) 
□ 



Proposition 3.10. Suppose that {^f}o<t<i is an ASP with generating law v under 
some measure P. Then 



dP 



■^t{Rt)-' 



(98) 



Tt 
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defines a probability measure Q for t e [0, 1). Furthermore, under Q {i,t}o<t<i is an 
n-dimensional gamma process such that 



n +_i 



Q[^t e dx] = n fj^^"'^ 

3.4 Independent gamma bridges representation 

The increments of an n-dimensional ASP are identical in law to a positive random 
variable multiplied by the Hadamard product of an n-dimensional Dirichlet random 
variable and a vector of n independent gamma bridges. For notational convenience, 
in this subsection we denote a gamma bridge defined over [0, 1] as {'y{t)} (instead of 
{7a}). 

For vectors x, y e M", we denote their Hadamard product as x o y. That is, 

xoy ^ (xiyi,...,XnynV- (99) 

Proposition 3.11. Given the value of $,g, the ASP process {^^} satisfies the following 
identity in law: 



{^t - is}s<t<i '= {R* D o 7j,<t<i, 



where 



1. D G [0, 1]" is a symmetric Dirichlet random variable with parameter vector (1 — 
s, . . . , 1 - sY; 

2. {7t} is a vector of n independent gamma bridges, each with activity parameter 
m — 1, starting at the value at time s, and terminating with unit value at time 

1; 

3. R* > is a random variable with law v* given by 

4- R* , D, and {7^} are mutually independent. 
Proof. Fix ki >1 and the partition 

s = 4<4<...<4 = i, (100) 

for i = 1, . . . , n. Define the non-overlapping increments {Aj^} by 

A,, = - ^ , (101) 
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for j = 1, . . . , /cj and i = 1, . . . ,n. The distribution of the vector 

A = (An,Ai2,..., Aifc,, 

A2I, A22, • • • , A2A;2, 



A„i,A„2,...,A„fcJ' (102) 

characterises the finite-dimensional distributions of the process {^j— ^s}s<i<i- It follows 
from the Kolmogorov extension theorem that the distribution of A characterises the 
law of {^i — ^<j}- Note that A are non-overlapping increments of the master GRB {Tt}. 
Thus, given A has a multivariate Liouville distribution with parameter vector 



i2 J.2 J.2 J.2 j2 j.2 

''I ~ ''0' '-2 ~ ''15 • • • 5 ''fc2 ~ ''A;2-l' 



^1 ^0'^2 ^L-i)"""' (103) 

and generating law 

i^*{A) = iystiA + Rs) (104) 
for t G (s, 1]. It follows from Fang et al. j^, theorem 6.9] that 

(A,i, . . . , A,,J^ '= R* AY„ for z = 1, . . . , n, (105) 

where (i) R* has law z/*, (ii) D = (Di,...,D„)^ has a Dirichlet distribution with 
parameter vector (1 — s, . . . , 1 — s)^, (iii) Yj G [0, 1]^' has a Dirichlet distribution with 
parameter vector {t\ — t^, . . . ,tl. — t^._i)^, (iv) Yi, . . . , Y„, R*, and D are mutually 
independent. 

Let {'j{t)}s<t<i be a gamma bridge with activity parameter m = 1 such that 7(5) = 
and 7(1) = 1. Then the increment vector 

(7(ti)-7(to),---,7(4J-7(tL-i)r (106) 

has a Dirichlet distribution with parameter vector (t\ — t^, . . . ,tl._ — tl._{)^ . Hence 
the increment vector fll06p is identical in law to Yj. From the Kolmogorov extension 
theorem, this identity characterises the law of {7(t)}. It follows that 

- ^f^s<t<i '= {R* Dat}s<t<u for ^ = 1, . . . , n, (107) 
which completes the proof. □ 
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3.5 Uniform process 

We construct a multivariate process from the ASP {^^} such that each one- dimensional 
marginal is uniformly distributed for every time t e (0, 1]. 

Fix a time t & (0, 1]. Each ^f'' is a scale-mixed beta random variable with survival 
function 

POO 

Ft{x)^ / {l-h/y[t,n-t])iy{dy) 

J X 

/•oo 

= / h_,/y[n-t,t]u{dy), (108) 

J X 

where Iz[q:, /3] is the regularized incomplete Beta function, defined as usual for z e [0, 1] 

by 



(109) 



The random variables 

Y^'^^F.i^i'^ t^l,...,n, (110) 

are then uniformly distributed. 

We now define a process {Y}o<t<i by 

Y,= (F,(e;')),...,F,(^(")))\ (111) 

By construction, each one-dimensional marginal Y^^^ is uniform for t > 0. 



4 Liouville process 

We generalise ASPs to a family of stochastic processes that we call Liouville pro- 
cesses. A Liouville process is a Markov process whose finite-dimensional distributions 
are multivariate Liouville. Liouville processes display a broader range of dynamics than 
ASPs. The one-dimensional marginal processes of a Liouville process are in general 
not identical. This generalisation comes at the expense of losing the direct connection 
to Archimedean copulas. 

Definition 4.1. Fix n e N+, n > 2, and the vector m e satisfying rui > 0, 
i — 1, . . . ,n. Define the strictly increasing sequence by 

Uq = 0, 

Ui = Ui_i + mi, fori = l,...,n. (112) 



26 



Then a process {$,t}o<t<i satisfying 



{it} 



law 



0<t<l 



t{u„—u„-i)+u„- 



u„-iJ ) o<t<l 



/or {ri}o<t<„„ a GRB with activity parameter m = 1, is an n- dimensional Liouville 
process. We say that the generating law of {Tt} is the generating law of {$^^\ and the 
activity parameter of {$,^} is m. 

Note that allowing the activity parameter of the master process to differ from unity 
in Definition l4.1l would not broaden the class of processes. Indeed, changing the activity 
parameter of the master process would be equivalent to multiplying the vector m by a 
scale factor. 

Let {^^} be a Liouville process with generating law v and parameter vector m. 
Each one-dimensional marginal process of {^^j is a GRB with activity parameter mj, 
and Definition 14. II ensures that is defined for t G [0, 1]. It is straightforward to adjust 
the definition so that a Liouville process is defined over an arbitrary closed interval. 

In the language of McNeil & Neslehova ll|, has a Liouville copula. That is, 
the survival copula of $,i is the survival copula of a multivariate Liouville distribution. 
We shall provide the transition law, moments and an independent gamma bridge rep- 
resentation of a Liouville process. We present the results as propositions. Proofs are 
omitted since they are similar to the proofs in Section 3. 

We define a family of unnormalised measures, indexed by t G [0, 1) and x G M+, as 



tpt{B]x) 



u{B), 

fT{l-t){z-x) 

B friz) 
r[T]e^ 



for 5 G -B(M) where T 



and 



r[T(i-t)] Js 

||m||. Again we write 

^>t{x) = V't([0,oo);x), 



(113) 



(114) 



(115) 



■t - \\<it 



(116) 



The process {Rt} is a GRB with activity parameter T. Given the law of Ri is 

V',(dr;i?3) 



z^si(dr) 



"^siRs 



117) 
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and law of Rt is 



exp( 



^.(i|x| 



T[T{t - s)] 



-dr 



(118) 



for t e {s, 1). 



Proposition 4.2. The Liouville process {^^} is a Markov process with the transition 
law given by 



rK(i - s)] 



anc? 



P[£te dy|C = x] 



MM 

^.fllxl 



n 

i=l 



_ T..'imi(t-s)-lp-(2/i-Xi) 



T[mi{t - s)] 



(120) 



where r(t) = 1 - (1 - t)/T, < s < t < 1, and B e B{R). 



Proposition 4.3. Fix s e [0, 1). T/ie /irst- and second-order moments of ^t, t e (s, 1], 

are 



E 



CO 



;2. Var 
3. Gov 



St ) St 



mi 
T 



mi{t - s) + 1 
T{t -s) + l 
mimj{t — s) 



T 



rrii 2 

fJ'2 



T{t - s) + 1 T{t- s) 



where 



fJ'2 = 



1-S 

{t-s){l + T{t-s)) 
{l-s){l + T{l-s)) 



E[{R^-R,y\R,]. 



The law of the increments of an n- dimensional Liouville process can be characterised 
by a positive random variable multiplied by the Hadamard product of an n-dimensional 
Dirichlet random variable and a vector of n independent gamma bridges. 

Proposition 4.4. Given the value of $,g, the Liouville process {^j} satisfies the follow- 
ing identity in law: 

{it - is}s<t<i ^= {R* D o 7t},<t<i, 



where 
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1. D e [0, 1]" has a Dirichlet distribution with parameter vector (1 — s)m; 

2. {'Jt} is a vector of n independent gamma bridges, such that the ith marginal 
process is a gamma bridge with activity parameter rrii, starting at the value at 
time s, and terminating with unit value at time 1; 

3. R* > is a random variable with law v* given by 
4- R* , D, and are mutually independent. 
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